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Transverse momentum spectra of electrons from Au+Au collisions at ^/s NN — 130 GeV have been 
measured at midrapidity by the PHENIX experiment at RHIC. The spectra show an excess above 
the background from photon conversions and light hadron decays. The electron signal is consistent 
with that expected from semi-leptonic decays of charm. The yield of the electron signal dN s /dy 
for pr > 0.8 GeV/c is 0.025 ± 0.004(stat.)±0.010(sys.) in central collisions, and the corresponding 
charm cross section is 380 ± 60(stat.)±200(sys.) fib per binary nucleon-nucleon collision. 

PACS numbers: 25.75.Dw 



In this letter, we report the first measurement of sin- 
gle electron spectra, (e + + e~)/2, in Au+Au collisions 
a t y/ s NN = 130 GeV at the Relativistic Heavy Ion Col- 
lider (RHIC). The measurement of single leptons at high 
transverse momentum (p? > 1 GeV/c) is a useful way 
to study heavy-quark production, an important probe of 
hot and dense matter created in high energy heavy ion 
collisions. Charm production is sensitive to the initial 
state gluon density Jl|,§. Nuclear and medium effects, 
such as shadowing and charm quark energy loss 
can be studied by comparison of charm production in 
AA, pA, and pp collisions. Measurement of charm is im- 



portant for understanding J/ip suppression (a proposed 
signal of the deconfinement phase transition |^^]) and 
the dilepton mass distribution in 1 < M;+;- < 3 GeV, 
where lepton pairs from charm make significant contribu- 
tions 0. Inpp collisions at the ISR (^/s = 30 — 63 GeV), 
production of single electrons was observed (e/ir ~ 10~ 4 ) 
for pt > 1 GeV/c [p|-pd|, and interpreted as evidence of 
open charm production In pp collisions at RHIC 

energies, the signal level is expected to be higher, since 
charm production increases with ^/s NJV faster than pion 
production. We recently observed suppression of high px 
pion production in Au+Au collisions at RHIC relative 
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to binary nucleon-nucleon (NN) collision scaling |L3| . If 
charm production scales with NN collisions, as expected 
in the absence of nuclear effects, the e/ir ratio will be 
even higher in Au+Au collisions at RHIC. 

Data used for this analysis were recorded by the 
PHENIX west-arm spectrometer @ (A(f> = 90° in az- 
imuth, \r]\ < 0.35 in pseudo-rapidity), which consisted of 
a drift chamber (DC), a layer of pad chambers (PCI), 
a ring imaging Cerenkov detector (RICH), and a lead- 
scintillator electromagnetic calorimeter (EMCAL). The 
trigger was provided by beam-beam counters (BBC) and 
zero-degree calorimeters (ZDC). ZDC and BBC signals 
were combined to select centrality: central (0-10%), pe- 
ripheral (60-80%), and minimum bias (0-92%) §J]. 

The analysis uses 1.23 M minimum bias events with 
vertex position \z\ < 30 cm. Charged particle tracks are 
reconstructed by the DC and the PCI with a momen- 
tum resolution Sp/p ~ 0.6% ffi 3.6% p (GeV/c). Tracks 
are confirmed by a matching hit in the EMCAL, which 
measures the energy E deposited with a resolution of 
8.2%/^S(GeV) © 1.9% for test beam electrons. Elec- 
tron identification is performed using the RICH and the 
EMCAL Q. The RICH is filled with 1 atm C0 2 and de- 
tects on average 10.8 photo-electrons per electron track, 
while a pion with p < 4.7 GeV/c produces no signal. It is 
required that at least 3 RICH hits are associated with the 
track and that their hit pattern is consistent with that of 
an electron track. After these cuts, a clear electron signal 
is observed as a narrow peak centered at E/p — 1.0. We 
select tracks in the peak as electron candidates. The E/p 
cut reduces hadron background, and removes conversion 
electrons created far from the vertex. A hadron deposits 
only a fraction of its energy in the EMCAL, and the 
momentum of an off-vertex conversion electron is recon- 
structed incorrectly. The remaining background, about 
10% of the electron candidates, is caused by accidental 
association of RICH hits with hadron tracks. The back- 
ground level is measured statistically by an event mixing 
method, and is subtracted from the yield. 

The electron acceptance (~ 7.4 % of dN/dy) and effi- 
ciency (~ 60 %) are determined using a detailed GEANT 
|l6f simulation, which satisfactorily reproduces the de- 
tector response. Additionally, a multiplicity dependent 
efficiency loss due to detector occupancy is evaluated by 
embedding simulated electrons into real events. This ef- 
ficiency loss is 27 ± 4% (4 ± 2%) for central (peripheral) 
collisions and has no significant px dependence. 

Fig. [l] shows the pt distributions of electrons in 
PHENIX for central, minimum bias, and peripheral col- 
lisions. Errors in the figure are statistical. The overall 
systematic uncertainty, which is the quadratic sum of sev- 
eral few percent effects, is about 11%. Expected sources 
of electrons are (1) Dalitz decays of tt° , r], rf , cj, and <j>, (2) 
di-electron decays of p, u, and <fi, (3) photon conversions, 
(4) kaon decays (if 0,± — > nev), (5) semi-leptonic decay 
of charm, and (6) other contributions such as bottom 



decays and thermal di-leptons. In this analysis, sources 
(l)-(4) are considered to be background. 

We have calculated the contributions from Dalitz 
and di-electron decays with a hadron decay generator. 
PHENIX has measured the pr distributions of ^ in 
0.2 < p T < 2.2 GeV/c @ and of tt° in 1 < p T < 
4 GeV/c |0. Since the ^ and w° data are consistent in 
the overlapping region, we fit a power law function to the 
combined data sets to determine the input 7r° spectrum 
for the decay generator. The px distribution of any other 
hadron h is obtained from the ir° spectrum by replacing 

Pt with Jp T + rrifr — rn^ - The shapes of the result- 
ing pt spectra of , p, and p agree with the PHENIX 
measurements ]l7j ] within 20%. In this parameterization 
h/ir° ratios approach constants at high pr- We assume 
the following asymptotic ratios to fix the relative normal- 
izations: 7]/tt° = 0.55, r/'/ir = 0.25, p/n° = lu/it = 1.0, 
<j>/TT° = 0.40. Except for the (f>, these ratios are taken 
from proton beam data of SPS, FNAL, and ISR experi- 
ments p8| , |l9| ]. The tj/tt ratio is consistent with a mea- 
surement in Pb+Pb collisions at SPS @. The 0/tt° ratio 
is based on the integrated ratio 4>/h~ ~ 0.02 in Au+Au 
collisions at ^Js NN — 130 GeV pi] ]. We assign to each 
ratio a conservative systematic uncertainty of 50%. 

Photon conversions are evaluated using a combination 
of the GEANT simulation and the hadron decay gener- 
ator. Since pt spectra of externally converted electrons 
are similar to those from Dalitz decay, the conversion 
spectra can be approximated by scaling the Dalitz de- 
cay spectra by an experiment specific factor, R con v = 
Conversion/Dalitz. R con v is evaluated using the GEANT 
simulation and is cross-checked by comparing the relative 
yield of reconstructed Dalitz and conversion pairs in the 
simulation and in data. The simulation shows that R CO nv 
has only a weak pt dependence, primarily due to the en- 
ergy dependence of the pair creation cross section. R CO nv 
is parameterized as (1.9 ± 0.2) x (1 - 0.0718 x p^ 76 ). 

Background from kaon decays is also evaluated using 
the GEANT simulation, and is found to be negligible. 

The upper panel of Fig. || shows the ratio of the mea- 
sured electrons to the calculated background versus pt 
for minimum bias events. The shaded region is the 
quadratic sum of systematic errors in the electron mea- 
surement and in the background. The latter includes 
uncertainties in the normalization and the shape of the 
7T° spectrum, in the h/n° ratios, and in R CO nv- A signif- 
icant electron excess above the background is observed 
for pt > 0.6 GeV/c. Central collisions show a similar ex- 
cess. The peripheral collision data sample lacks sufficient 
statistics to reveal a signal in this analysis. 

Fractional contributions to the background are shown 
in the lower panel of Fig. pi More than 80% of the back- 
ground is from ir° decay, directly from the Dalitz decay 
or indirectly from photon conversion. The ir° spectrum is 
well constrained by the PHENIX measurement. The next 
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most important background source is r\ decay. Given the 
assigned systematic error, the upper limit of the high px 
asymptotic 77/71"° ratio is 0.83. Since this ratio, corrected 
for feed-down, would imply that the primary ij/tt ~ 1, 
this provides a conservative limit on contributions from 
77's. Contributions from all other hadrons combined are 
only a few percent of the total. 

Background-subtracted electron spectra are shown in 
Fig. ||. The error bars on the data points represent the 
statistical errors, while the systematic error due to the 
background subtraction is indicated by brackets. The 
integrated yield of the electron signal dN e /dy for px > 
0.8 GeV/c is 0.025 ± 0.004(stat.)±0.010(sys.) for central 
collisions and is 0.0079 ± 0.0006(stat.)±0.0034(sys.) for 
minimum bias collisions. 

Semi-leptonic decay of charmed hadrons is an expected 
source of the electron signal. We use the event generator 
PYTHIA to estimate electron spectra from charm de- 
cay. We tuned the parameters [^3| of PYTHIA such that 
charm production data at SPS and FNAL j24j] and single 
electron data at the ISR [§-11 are well reproduced. The 
charm production cross section in pp collisions from this 
PYTHIA calculation is a c - c = 330/jb at V« =130 GeV. 
The electron spectrum in Au+Au collisions is then calcu- 
lated as EdN e /dp 3 = Taa x Eda e /dp 3 , where Eda e /dp 3 
is the electron spectrum from charm decay calculated 
with PYTHIA, and T AA (listed in Table @) is the nuclear 
overlap integral calculated from a Glauber model Jl3| . 
The calculated electron spectra shown in Fig. [5] are in 
reasonable agreement with the data. 

Before attributing the entire electron signal to open 
charm decays, it is necessary to quantify contributions 
from other possible sources. An analogous PYTHIA esti- 
mate of the bottom decay contribution is shown in Fig. |^. 
It becomes significant only above the measured pr range. 
Expected contributions from J/\l/ and Drell-Yan are neg- 
ligible. In Pb+Pb collisions at SPS, direct photons [|2y] 
and an enhanced yield of low mass di-leptons [^5]| have 
been reported. If these are due to thermal radiation 
from hot matter, an even larger production is expected at 
RHIC energies, and can contribute to the electron signal. 
Since p — ► e + e~ contributes less than 1% to the calcu- 
lated background as shown in Fig. ||, and since the dom- 
inant source of thermal di-leptons is 7r + n — > p — > e + e~ 
p6| , a significant contribution from thermal di-leptons is 
unlikely. There are several predictions for direct photons 
at RHIC energies The conversion electron spec- 

trum calculated from a prediction in ref. |2^] is shown in 
Fig. U for central collisions. It could explain 10-20% of 
the signal, with large theoretical uncertainties. 

Neglecting these other possible sources and assuming 
that all the electron signal is from charm, we derive the 
charm cross section corresponding to the electron data. 
We fit the charm electron spectrum from PYTHIA to the 
data for pt > 0.8 GeV/c, and obtain the rapidity density 
dN c5 /dy\ y= o and the total yield N ce of open charm. They 



are then converted to cross sections per NN collision: 
da cS /dy = (dN c£ /dy)/T A A and er cg = N cE /T A a- Results 
are shown in Table Q. The systematic error is a quadratic 
sum of many sources. For central collisions, they are 
background subtraction (±44%), uncertainties in the 
PYTHIA calculation (±11% from < k T >= 1.5 ± 0.5, 
±13% from D+/D = 0.65 ± 0.35, ±8% from PDFs), 
fit range (±18%), and Taa (±7%). Note that any finite 
contribution from neglected sources would reduce the de- 
rived charm cross section. Without nuclear or medium 
effects in charm production, a cS per NN collision should 
be independent of centrality. Within uncertainties, our 
data are consistent with this expectation, in possible con- 
trast to the attribution of increased charm production as 
the source of enhanced di-muon production reported in 
Pb+Pb collisions at SPS f§. 

The single electron signal yield (divided by Taa to give 
the cross section per NN collision) and the derived charm 
cross section are compared with single electron data of 
ISR experiments and charm data of fixed target experi- 
ments 23] in Fig. [|. Cross section curves calculated with 
PYTHIA, which has been tuned to the charm data and 
the ISR electron data, and a charm cross section curve 
from a next-to-leading order (NLO) pQCD calculation 
p(| are also shown in the figure. Our data are consistent 
with both of the calculations within large uncertainties. 

In conclusion, we have observed single electrons above 
the expected background from decays of light hadrons 
and photon conversion in Au±Au collisions at y/s NN = 
130 GeV. The observed signal is consistent with semi- 
leptonic decay of charm. The forthcoming high statistics 
Au±Au data and pp comparison data at full RHIC en- 
ergy (V s n« — 200 GeV) will be useful to clarify the na- 
ture of the single electron signal and to better determine 
heavy-quark production in Au±Au collisions at RHIC. 
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Au+Au -> (e+e")/2 + X (\|s^,=130GeV) 
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FIG. 1. Transverse momentum spectra of electrons in 
PHENIX from Au+Au collisions at ./s7^7=130 GeV. 
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FIG. 2. Ratio of the electron data to the calculated back- 
ground as a function of pr in minimum bias collisions (upper 
panel) and relative contributions to the background from var- 
ious sources (lower panel). The curves for u) and <j) show the 
sum of the Dalitz and the di-electron decay modes. 
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FIG. 3. The background-subtracted electron spectra for 
minimum bias (0-92%) and central (0-10%) collisions com- 
pared with the expected contributions from open charm de- 
cays. Also shown, for central collisions only, are the expected 
contribution from bottom decays (dashed) and the conversion 
electron spectrum from a direct photon prediction (dotted). 
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FIG. 4. Single electron cross sections da e /dy\ y= o of this 
measurement and ISR experiments [p|, pj]j3l| ] are displayed 
(bottom of fig., right-hand scale) with charm decay contribu- 
tions calculated with PYTHIA. Open and filled symbols are 
for 1.0 < pr < 1.4 GeV/c and pr > 1.4 GeV/c, respectively. 
The derived charm cross section of this measurement is com- 
pared with charm cross sections from SPS/FNAL experiments 
(top of fig., left-hand scale). The thick curve and the shaded 
band represent the charm cross section in the PYTHIA model 
and in a NLO pQCD calculation |30|, respectively. 



TABLE I. Charm cross section per NN collision derived 
from the single electron data for central (0-10%) and mini- 
mum bias (0-92%) collisions. The first and second errors are 



statistical and systematic, 


respectively. 




Centrality T^mb -1 ) 


dcr C c/dy\y=o(txh) 


<7 C c(/ib) 


0-10% 22.6 ± 1.6(sys.) 
0-92% 6.2 ± 0.4(sys.) 


97 ± 13 ±49 
107±8±63 


380 ± 60 ± 200 
420 ± 33 ± 250 
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